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Abstract 

We discuss the electromagnetic energy-momentum distribution and the me- 
chanical forces of the electromagnetic field in material media. There is a 
long-standing controversy on these notions. The Minkowski and the Abra- 
ham energy-momentum tensors are the most well-known ones. We propose a 
solution of this problem which appears to be natural and self-consistent from 
both a theoretical and an experimental point of view. 
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I. INTRODUCTION 



The discussion of the energy-momentum tensor in macroscopic electrodynamics is quite 
old. The beginning of this dispute goes back to Minkowski [1], Abraham [2], and Einstein 
and Laub [3]. Good reviews of the problem can be found in [4-7], to mention but a few. 
Nevertheless, up to now the question was not settled and there is an on-going exchange of 
conflicting opinions concerning the validity of the Minkowski versus the Abraham energy- 
momentum tensor, see [8] for a recent discussion. Even experiments were not quite able to 
make a definite and decisive choice of electromagnetic energy and momentum in material 
media. 

Here we propose the solution of the problem. 

Our basic notations and conventions are as follows. We are using international SI units 
throughout. Correspondingly, e , /i are the electric and the magnetic constant (earlier 
called vacuum permittivity and vacuum permeability). The Minkowski metric is = 
diag(c 2 , — 1, — 1, — 1). Latin indices from the middle of the alphabet label the spacetime 
components, k, . . . = 0, 1, 2, 3, whereas those from the beginning of the alphabet refer to 
3-space: a, b, c, . . . = 1, 2, 3. 



Our solution is as follows. The new electromagnetic energy-momentum in an arbitrary 
medium reads: 



The electromagnetic field strength = (E, B) is composed of the electric and magnetic 
3-vector fields. Componentwise, Eq.(l) describes the energy density of the field 



II. NEW ENERGY-MOMENTUM TENSOR 
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its energy flux density (or Poynting vector) 
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T a = s a = — [E x Bf , (3) 

its field momentum density 

T a ° = - Pa = e [B x E] , (4) 

and the Maxwell stress tensor 

As one can immediately notice that (1) has the same form as the vacuum energy- 
momentum tensor. However, its physical content is different, which follows from the fact 
that Fij = (E, B) satisfies the macroscopic Maxwell equations in matter: 

V x H - D = j free , V • D = p free , (6) 
VxE + B = 0, V-B = 0. (7) 

Here p free and j free are the densities of the free (or, in other terminology, the "true" or "ex- 
ternal") charge and the free current. The fields D and H represent the electric and magnetic 
excitations (other names are "electric displacement" and "magnetic field intensity"). The 
field strengths (E, B) and the excitations (D, H) are related by means of the equations 

e E = D-P, and — B = H + M. (8) 

The polarization P and the magnetization M represent the bound (or "polarizational" ) 
current and charge densities inside the material medium: 

j mat = P + V x M, p mat = - V • P. (9) 

The density of the mechanical (ponderomotive) force acting on matter is determined as 
the divergence of the energy-momentum tensor 

fi = djTJ. (10) 



Differentiating (1) and using the Maxwell equations (6), (7) and the equations (8) and (9), 
yields: 



/o = « + V-s=-E-j tot , 

fa = ~ Pa + d b S a b = P t0t E a + [j t0t X B 



(11) 

(12) 



Here the total charge and current density are p tot = p frcc + p mat and j tot = j frcc + j mat . 

This result is quite natural and physically clear. The electromagnetic field affects matter 
by means of the two Lorentz forces (12): One acts on the free charge and current (p free , j free ) 
(on the conductive current, for example), and another force acts on the bound charge and 
current (9). The latter have also a direct physical meaning in the macroscopic (Lorentz type 
averaging) framework and in the microscopic approaches, see Hirst [9], e.g. The temporally 
and spatially varying polarization and magnetization give rise to the electric and magnetic 
fields, like the free charges and currents do. Conversely, the bound charges and currents 
should also feel the electromagnetic field in the same way as the free charges and currents 
do. 

The representation of the total electromagnetic force as the sum of the two terms with the 
clear-cut physical content (12) suggests a natural step to split the original energy-momentum 
(1) into the corresponding sum of the two energy-momentum tensors which are associated 
with the free and bound charge/current, respectively. Using (8) in (1), we find 



T>.3 — trpj I brpj 



(13) 



where we introduce the free-charge energy- momentum and the bound-charge energy- 
momentum tensors as 



( 



u 



{ Pa^ 



fc b 









~ h Pa 




bob 



(14) 



The components of the free-charge energy-momentum read explicitly 

l u = - (E • D + B • H) , 



(15) 



f s = E x H, f p = D x B, 



{ S a b = E a D b + B a H b - 1 5* (E • D + B • H) , 

whereas the components of the bound-charge energy-momentum are 

h u = -(-E-P + B-M), 

b s = ExM, b p = -PxB, 
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h S a b = - E a P b + B a M b - ^ 6* (- E • P + B • M) . 



(16) 
(17) 

(18) 
(19) 
(20) 



One straightforwardly recognizes the tensor with the components (15)-(17) as the 
well-known Minkowski energy-momentum tensor. 

Similarly to (10), the divergences f /« — d/T^ and b fi = d^Tj determine the force 
densities. By construction, the total 4-force density is the sum /j = f /j + h f\. Explicitly, we 
have for the 3-force densities 



i f a = p hcc E a + [j 



free x g 



h f a = p™ t E a + j 



mat x g 



-x n 



+ x a . 



Here, in the Cartesian coordinates, 

X a = \ (E b d a D h - D b d a E b + H b d a B b - B b d a H b ) . 



(21) 
(22) 

(23) 



In particular, X a — \ e E b E c d a e bc + p: H b H c d a p bc for the linear material laws D a = 
Eo e ab E b and B a = /x p ab B b . This extra term vanishes for homogeneous media. 

Thus indeed, tensors in the sum (13) are associated with the two different types of 
charges and currents in the material medium. The divergences produce, essentially, the 
two independent Lorentz forces acting separately on the free and on the bound charge and 
current. 



III. SOME PROPERTIES AND APPLICATIONS 



After all these preliminaries and formal derivations, we are in a position to discuss 
the physical properties of the energy-momentum (1). At first, some remarks about (1) in 
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comparison to the Minkowski and Abraham tensors. Many authors (see the discussions in 
[5-7]) pointed to a clearly unphysical result produced by the Minkowski energy- momentum: 
in the absence of free charges and currents, a homogeneous medium appears to be always 
subject to the zero electromanetic force. This fact was usually taken in favor of the Abraham 
tensor which predicts an extra, so called Abraham force. However, the energy-momentum 
(1) does not suffer from such a deficiency. Even when the free charge and current densities 
are vanishing, the total force is, in general, non-trivial in view of the presence of polarization 
charge and current. Moreover, as compared to the rather ad hoc choice of the Abraham 
force, the mechanical action on the bound charge and current is in all cases described, using 
(1), by the well-known Lorentz force (12). 

Furthermore, the Minkowski tensor is asymmetric which is obvious from the comparison 
of the energy flux f s and the field momentum f p in (16). Usually, this fact was also taken in 
favor of the Abraham tensor, which is symmetric. At the same time, despite its symmetry, 
the structure of the Abraham tensor is defined in a rather ad hoc manner with opaque 
physical motivations. In contrast to this, the energy-momentum (1) is naturally symmetric 
and the electromagnetic field momentum (4) is related to the energy flux (3) as p = s/c 2 . 
This is the famous Planck relation which generalizes the Einsteinian Am = AE/c 2 relation 
to field theory. The interesting discussion of von Laue [10] was concentrated mainly around 
this point. 

As a simple application, let us consider the propagation of an electromagnetic plane 
wave from the vacuum into a dielectric medium with /j, — 1 and refractive index n = ^fe. 
More exactly, like in the previous discussions [11,6,7], we will confine ourselves to the case 
of a gaseous medium consisting of heavy atoms. We assume normal incidence on a plane 
boundary and we recall the reflection and transmission coefficients R = (n — l)/{n + 1) and 
T = 2/(n + 1), respectively. Then, for incident and reflected waves in vacuum, we find for 
the mean field energy (15) and momentum (16), if averaged over one period, 

u = e ^(l + R 2 )=e \E f±^ 2 , (24) 
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c 2 _V ^ " y " c '~ UI {l + nf 



- £ ° |Eo1 :i- J R 2 )k = ^|E | 2 -^-.k. (25) 



On the other hand, within the dielectric, for the transmitted wave, Eqs. (15), (18) and (16), 
(19) yield 

l + n 2 \E \ 2 ^ 2 1 + n 2 

M = £ ° — 2~ T =£olE ° l (1+^)2' (26) 

P = ^^T 2 k=^W-^-k. (27) 
c 2 c (l + n) 2 

Here i?o is the amplitude of the electric field and k the unit wave vector which specifies the 
direction of propagation. Comparing the above formulas, we see that both, the total energy 
and the total momentum calculated on the basis of (1), are conserved on the passage of the 
wave into the medium. This conclusion is in a complete agreement with the previous studies 
[11,6,7]. Moreover, it can be supplemented by a far more detailed analysis of wave propaga- 
tion in a gaseous media that has demonstrated [11,6,7] the plausibility of the ("Abraham") 
field momentum Ex H/c. However, since all these studies were confined to dielectrics with 
\i = 1, the arguments presented in [11,6,7] in actual fact give support to the field momentum 
(4) likewise. 

IV. WALKER AND WALKER EXPERIMENT 

Finally, let us discuss the direct experimental confirmations of the energy-momentum 
(1). For this purpose, as a first example, we recall the measurements of Walker & Walker 
[12] of the force acting on a dielectric disk placed in crossed oscillating electric and magnetic 
fields, see Fig. 1(a). The scheme of the experiment is as follows: The symmetry of the 
problem suggests to use cylindrical coordinates (r, (p,z). A small cylinder is made of barium 
titanate with e = 3340 and fj, — 1. Its height in z-direction is I ~ 2 cm, internal radius 
r\ ~ 0.4 cm, and external radius r 2 ps 2.6 cm. This disk is suspended between the poles of 
an electromagnet which creates the harmonically oscillating axial magnetic field. Besides, 
the oscillating field, with the phase difference of tt/2, a radial electric field is created by 
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means of an alternating voltage applied between the inner and the outer cylindrical surfaces 
of the disk. The oscillation frequency of the fields is rather low, namely, u = 60 Hz. As a 
consequence, everywhere in the disk rujn/c ~ zum/c ~ 10~ 7 <C 1. Walker & Walker [12] 
measured the torque around the z-axis produced by the electromagnetic force. 

Let us derive the theoretical value of the torque which is given by the volume integral 

N z = f dv[rxi] z = ( dvrf^. (28) 

J disk J disk 

Since there are no free charges and currents in the dielectric, the Lorentz force (12) reduces 
only to the second term determined by the bound charges and currents. One can check 
that the Maxwell equations (6), (7), together with the constitutive relations D = ee E and 
B = /ioH, are solved by the electric and magnetic field configuration 

E = sin( ^H^^ e ' + T^)- (29) 

B = cosM) (b e, - rlog ^ /ri )c2 (30) 

These approximate formulas are valid with very high accuracy due to the fact that 
rujn/c ~ zum/c ~ 10~ 7 <C 1 everywhere in the disk. As usual, (e r ,e v ,e z ) denote the local 
orthonormal frame vectors of the cylindrical coordinate system. Here Bq is the magnitude 
of the oscillating axial magnetic field and Uo is the amplitude of the voltage applied between 
the inner (r = r-i) and the outer (r = r 2 ) cylindrical surfaces of the disk, AU = U sm(ujt). 

The bound charge and current densities (9) are straightforwardly found to be p mat = 
and j mat = P = Eq(£— 1)E. Correspondingly, the Lorentz force turns out to be f = j mat x B = 
e (e - 1)E x B. Substituting (29) and (30) in (28), we obtain the torque 

N z = -eo(e - l)nl(rj - r?) = U f Q " , cos 2 ^). (31) 

log(r 2 /ri) 

This result was experimentally confirmed by Walker & Walker [12]. Although the authors 
of [12] apparently noticed that the torque measured fits the understanding of the electro- 
magnetic force as the Lorentz force for the polarization current (in agreement with our 
approach), they ultimately claimed that their experiment confirms the Abraham force. The 
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theoretical explanation presented in [12] was based on the idea that the Maxwell stress S a b 
caused the unusual surface drag of the disk. However, when computing the force, see (12), 
the contribution of the stress dbS a b should be complemented by the momentum term — p a ; 
the latter is missing in [12]. The alternative explanation [5] is based on the computation of 
forces on the metal coatings of the cylinder. This yields the result where the factor (e — 1) 
above is replaced by e. Since the dielectric matter in question has e = 3340, this experiment 
thus cannot be treated as the critical test for the different energy-momenta. 



As another example we consider the experiment of James [13] which is in many respects 
very similar to the one of Walker & Walker. James, see Fig. 1(b), had also placed a disk 
into crossed electric and magnetic fields. The small cylinders were made of a composition 
of nickel-zinc ferrite with /j, — 16 or 43 and e ~ 7. Like in [12], the radial electric field 
was created by means of an oscillating voltage applied between the inner and the outer 
cylindrical surface of the disk. However, instead of an axial magnetic field, an azimuthal 
magnetic field was produced inside matter by an alternating electric current in a conducting 
wire placed along the axis of the disk. The resulting field configuration reads: 



These formulas hold true in the approximation zu u n/c ~ zuiii/c ~ 1CT 5 <C 1 which is 
fulfilled in James' experiment everywhere in the cylinders with r±, r 2 and the length / of order 
1-3 centimeter. Here Jo is the amplitude of the alternating current I — Iq sin(cjjt) along the z- 
axis, whereas Uq gives, as before, the amplitude of the oscillating voltage, AU = Uq sm(uj u t). 
The frequencies u>i and u u are different and are varied in the course of the experiment between 
10 and 30 kHz. Since this experiment, unlike [12], covers also magnetic media, we display 
the nontrivial permeability //. One can check that (32), (33) satisfy the Maxwell equations 



V. EXPERIMENT OF JAMES 




(32) 



(33) 
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(6), (7). 



A. Electromagnetic force in James' experiment 

James [13] measured the force T z acting along the axis of the disk in the crossed fields 
(32), (33). Let us derive the theoretical value of this force by using the general expression 
for the force density (12). There are no free charges and currents inside matter, p frec = 
and j frcc = 0. By substituting (32), (33) into (9), we find p mat = and 

•mat / ( U) U U Q C0S{UJ) /i/i / . , .s 2 \ f 0A \ 

J = - 1) [ rlog(ra/ri) + M^t) ^ z) e r - (34) 

The total force is obtained as the integral of the force density 



;mat 



x B 



(e/x - l)e r nno U I 



■ oj u sin(c^jt) cos(w M t) 



-2tt log(r 2 /ri) 

(nf/ /clog(r 2 /ri)) 2 u 2 u cos 2 (u u t) z + (fi/i I /27r) 2 u 2 sin 2 (cu;t) z e z (35) 



over the volume of the disk: 

T z =\ dv [j mat x B] 2 = ~ Uu sin (^ t ) C os(^). (36) 

J disk CT 

According to James [13], we choose uJi = u u ±u , with u the mechanical resonance frequency 
of the disk, and find the final expression for the force 

2c 2 Uu sm ^ ^- ( 37 ) 

This theoretical prediction was actually verified in the experiment of James [13]. 



B. Minkowski and Abraham forces in James' experiment 

In the isotropic case under consideration, e ab = 5 ab e and fi ab = 5 ab \x. With the 
free charges and currents absent, the Minkowski 3- force density (21) reduces to the last 
term which contributes only at the ends of the cylinder. Since the permittivity has 
the constant value of e ^ 1 inside the body, i.e. for — 1/2 < z < 1/2, and drops to 
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e = l outside of that interval, the derivative of such a stepwise function reads d z e(z) = 
(e — 1) [8(z + 1/2) — 5(z — 1/2)}. Similar relation holds for the derivative of the permeability 
function fx(z). Correspondingly, we find for the Minkowski force 



/ 



Tf = - I dvX z 



disk 

T2 



1 

2 J disk 



dv 



e Q E 2 d z e + fi H 2 d z pL 



= 7T £ 2 drr{(e - l)e [E 2 (l/2) - E 2 (-l/2)] + (/i - l)/i [H 2 (l/2) - H 2 (-l/2)] } . (38) 

Substituting the squares of the electric and magnetic fields (32) and (33) [note that H = 
B//i/i ], we obtain: 



T. 



M 



IU I 



[li(s — 1) uJi sin(c<j u t) cos(cjjt) + e(fx — 1) u u sm(ujit) cos(u> u t)] . 



(39) 



The Abraham 3-force density differs from the Minkowski expression by the so-called 
Abraham term [see [5], eq. (1.6) on page 140, for example]: 



f A = f M + ^|(ExH). 

<r at 

It is straightforward to evaluate the last term. Using (32) and (33), we get 



(40) 



^ExH = 4 

c z fir z 



fifxo U I 



2ir log(r 2 /ri) 



sin(u u t) sin(c<Jjt) — ( / u / u /o/2vr) 2 UJi sin(cjjt) cos(cjjt) z 



- (nU /c\og(r 2 /ri)) 2 lu u sin(cj u t) cos(cj u t) z 



e 2 . 



(41) 



Taking the time derivative and integrating over the body, we find an additional contribution 
to the total force: 



, ue — 1 d 
dv P - -(Ex H), 



disk C 2 dt 

(fie - 1) IUqIq 



[uJi sin(uj u t) cos(cjjt) + u u sin(wjt) cos(o> u t)] . 



(42) 



In James' experiment, we put ut = uj u ± u and select only the component varying with 
the mechanical resonance frequency of the body, uoq. Then (39) and (42) yield the Minkowski 
and the Abraham forces: 

IU I 



2c 2 (v-£)u u ]sm(u Q t), 



2c 2 



[(1 - ji) ujq ^ (// - e) uj u ] sin(a;ot). 



(43) 
(44) 
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C. Theories versus experiment 



All the theoretical expressions for the electromagnetic force look similar: compare 
(37) with (43) and (44). However, the crucial difference is revealed when we recall that 
James measured not the force itself but a "reduced force" defined as the mean value 
| [JF z (c<j u , Ui = uj u + u ) + T z {oj u ,ijJi = 0J U — coo)], see eq. (9) on page 60 of James' thesis 
[13] and the footnote on page 158 of [5]. With high accuracy, James observed the vanishing 
of the reduced force in his experiment. This observation is in complete agreement with the 
theoretical derivation (37) based on our new energy- momentum tensor, whereas both, the 
expressions of Minkowski (43) and of Abraham (44), clearly contradict this experiment. 

The explanation proposed in [5] in support of the Abraham force appears to be incon- 
sistent mathematically and misleading physically. Namely, the computation of the force is 
reduced in [5] to the evaluation of the surface integral of the Maxwell stress in the vac- 
uum "just outside the disk". However, instead of the usual continuity of the tangential 
electric field, an unsubstantiated matching condition was introduced for E on the boundary 
(z = ±1/2) in order to find the fields outside the disk. Such a derivation [which yields a result 
different from (44) above] cannot be considered to be a satisfactory theoretical explanation. 

To begin with, there is not any good reason why one should replace a well-defined volume 
integral for the total force by a surface integral. Formally, this is allowed, of course, but as 
soon as we know the fields inside the body everywhere, see (32) and (33), we can proceed 
directly by constructing the explicit expressions of the force densities (21), (22), and (40) 
and then straightforwardly find the corresponding volume integrals. There is no logical need 
to perform an auxiliary computation in order to find the vacuum fields "just outside" the 
body, which appears to be a separate nontrivial problem. Provided the latter problem is 
solved correctly, we anticipate that the final result would agree with our (44). And certainly, 
one should use the standard matching conditions since this amounts to nothing else than to 
apply Maxwell's equations in a thin neighborhood near the surface. Accordingly, as to the 
matching of the electric field, we can only impose (as usual) the continuity of the tangential 
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components of electric field. Imposing a different discontinuity condition [as was done in eq. 
(3.17) of ref. [5], e.g.] is tantamount to assuming that Maxwell's equations are violated near 
and across the boundary. 

In our theoretical analysis, we used the field configuration (32), (33) which is valid 
inside of the cylinder. Near the ends, strictly speaking, one should take into account the 
deformation of the fields. However, it is well known that the corresponding corrections are 
confined to the regions very close to the ends. More exactly, the most important point is 
that the resulting end corrections for the total force are not proportional to the length of 
the cylinder. In other words, such end corrections (provided one computes them carefully) 
obviously would not compensate the reduced force of Minkowski (43) and of Abraham (44), 
which are both proportional to the length /. It is worthwhile to note that the end corrections 
were never taken into account in the previous analyses [13,5], and we use here precisely the 
same field configuration (32), (33) as in [13,5]. 

VI. DISCUSSION AND CONCLUSION 

Let us summarize our results. In the present paper we gave evidence that the correct 
energy-momentum tensor of the electromagnetic field in material media is described by (1). 
This tensor is symmetric and satisfies Planck's field-theoretical generalization p = s/c 2 
of Einstein's formula Am = AE/c 2 . The corresponding electromagnetic force turns out 
to be the Lorentz force acting on the free and bound charge and current densities. The 
energy-momentum (1) can be naturally represented as a sum (13) of the Minkowski energy- 
momentum and the bound-charge energy-momentum tensor. 

Our derivations here are in fact motivated by our axiomatic approach to classical elec- 
trodynamics [14] in which the Lorentz force represents one of the fundamental postulates of 
the scheme. In particular, if one starts from the Lorentz force equations (11) and (12) and 
reverses the order of the equations, one finally derives the energy-momentum tensor (1) that 
we first introduced without preliminary explanations. Besides the evidence of the general 
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validity of the Lorentz force axiom for point particles, a careful analysis of the wave prop- 
agation in material media as well as a proper interpretation of the experiments by Walker 
& Walker and by James, give further support to this basic cornerstone of classical elec- 
trodynamics. In our discussion we did not touch the electro- and magnetostriction effects 
because their consideration requires a more detailed specification of the internal mechanical 
properties of the medium. Moreover, in most cases the overall electro- and magnetostriction 
effects are balanced and are not directly observable. 

At the present level of understanding, we can thus conclude that the tensor (1) passes 
the theoretical and experimental tests and qualifies for a correct description of the energy- 
momentum properties of the electromagnetic field in macroscopic electrodynamics. 

As we have learned recently, the same energy-momentum tensor was introduced by 
P. Poincelot [15] who insisted on the equal physical treatment of the free and the polar- 
izational charges and currents. Such an equality is natural in our axiomatic approach to 
electrodynamics [14]. 
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FIG. 1. (a) Walker & Walker experiment; (b) James's experiment. 
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